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We report the experimental observation of conductance quantization in graphene nanoribbons, 
where ID transport subbands are formed due to the lateral quantum confinement. We show that 
this quantization in graphene nanoribbons can be observed at temperatures as high as 80 K and 
channel lengths as long as 1.7 /xm. The observed quantization is in agreement with that predicted 
by theoretical calculations. 



The isolation of graphene [T], a single sheet of graphite, 
has incited numerous studies because of fundamental 
physics interests |5] |3] and promising applications for 
carbon-based electronics . As a truly two-dimensional 
system and a zero-gap semiconductor where the elec- 
trons and holes behave as massless fermions, graphene 
possesses distinctly different transport properties from 
that of conventional 2D and 3D electronic materials. 
However, in order to utilize their remarkable electri- 
cal characteristics in certain applications, it would be 
highly desirable to produce a band gap in graphene, 
and therefore, intense efforts are being made to explore 
the properties of low-dimensional (ID and OD) graphene 
nanostructures [5J [SJ H 1^ . 

By patterning graphene into a narrow ribbon struc- 
ture, the carriers are laterally confined to form a quasi- 
one-dimensional (ID) system, similar to the case of car- 
bon nanotubes. Due to the linear dispersion relation 
E =^ Vfhk and high Fermi velocity (vf lO^m/s) in 
graphene, the quantization energy of graphene nanorib- 
bons (GNRs) can be substantially larger than that of 
conventional semiconducting materials of the same di- 
mension and parabolic dispersion. The formation of ID 
subbands in GNRs is expected to yield an energy gap for 
certain ribbon widths and crystallographic directions [H]. 
In terms of transport, this quantum confinement can 
also lead to quantized conductance, which is one of the 
most important transport characteristics of mesoscopic 
physics. 

Recently, experimental studies on GNRs have revealed 
a thermally-activated conductivity [3 [6] , suggesting the 
presence of a size-dependent energy gap. However, there 
remains some controversy in understanding the transport 
behavior in realistic GNR devices, where the imperfect 
and unknown edge configurations may lead to Columbic- 
blockade type transport at low temperatures [8] . Fur- 
thermore, an energy gap could also be induced by sub- 
strate interactions |10j. Despite numerous theoretical pre- 
dictions on GNRs [TTl [T^ ITB] , conductance quantization 
has yet to be reported experimentally. 

In this Letter, we present electrical transport measure- 
ments for GNR devices with lateral widths of 30 nm, and 
report conductance plateau features in GNRs at tem- 
peratures as high as 80 K. By modulating the Fermi en- 
ergy by a back gate, we observe conductance quantiza- 
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FIG. 1: (a) Schematics of the GNR device, where the 
graphene is contacted by two Pd metal leads, and the p-doped 
Si substrate acts as the gate electrode, (b) AFM image of 
GNR devices with different channel lengths, (c) Resistance 
of a GNR device measured as a function of gate voltage before 
and after the HSQ layer is removed, showing the impact of 
HSQ on electrical behaviors. Both measurements were per- 
formed in vacuum after annealing at 135°C. (d) Conductance 
of a GNR device {W = 30 nm and L = 850 nm) measured as 
a function of gate voltage at 300 K and 80 K. The bias voltage 
is 10 mV. 

tion plateaus for electron and holes in the same device, 
both possessing comparable transmission probability for 
each ID conduction mode. Our results provide the di- 
rect experimental evidence of quantum size confinement 
effects and the formation of subbands for ID graphene 
nanostructures . 

The GNR devices in our study, schematically shown 
in Fig.[l|a), were fabricated from mechanically-exfoliated 
graphene sheets on a p-doped Si substrate covered with 
300-nm thick Si02. Raman spectroscopy and atomic 
force microscope (AFM) measurements were employed 
to identify single-layer graphene. For the details of the 
device fabrication, we refer to Ref. [3 . The source 
and drain contacts are made of Pd, and the Si substrate 
acts as the back gate. Narrow ribbons of graphene are 
formed by oxygen plasma RIE etching using a patterned 
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HSQ (hydrogen silsesquioxane) layer as the protective 
mask. Fig.[ljb) shows the AFM image of an array of 
30-nm-wide GNRs after this HSQ layer is removed in hy- 
drofluoric acid (HF) solution. We found that electrical 
properties of the GNRs are strongly affected by this HSQ 
layer, which needs to be removed in order to reveal their 
intrinsic transport properties and, more importantly, the 
conductance quantization behaviors at low temperatures, 
as discussed below. Fig.[l](c) shows the resistance as a 
function of gate voltage Vg of a GNR device before and 
after the HSQ layer is removed. We note that the pres- 
ence of HSQ leads to a significant positive shift of the 
Dirac point voltage and a reduction of carrier mobility 
in graphene. In contrast, after removing the HSQ layer, 
the GNR exhibits ambipolar behavior with the resistance 
maximum associated with the Dirac point occurring at 
~ V, indicating a negligible unintentional doping in 
the final device. 

Transport measurements of GNRs were performed in 
vacuum (~ 10~^torr), and the devices studied all pos- 
sess a Dirac voltage near Vg ~ OV after annealing at 
135°C. Fig.[l];d) shows the conductance G of a GNR de- 
vice (L = 850 nm and W = 30 nm) measured as a func- 
tion of gate voltage Vg at a dc bias of 10 mV. At room 
temperature, the conductance curve G{Vg) in Fig.jlj^d) 
resembles that of bulk graphene [5] , showing the char- 
acteristic "V" shape that reflects symmetric hole and 
electron transports at negative and positive gate volt- 
ages, respectively. While the overall conductance curve 
G{Vg) of the GNR displays little variation as tempera- 
ture decreases, several plateau features start to appear 
in the measured G{Vg) curve and become apparent for 
T < 100 K (see Fig.jljd)). The shght asymmetry in the 
slope of the n and p branches is likely associated with the 
gate oxide hysteresis. Nevertheless, we note that these 
conductance plateaus are observed in both electron and 
hole branches with nearly the same conductance values 
and an equal spacing (see dashed line in Fig.[T]jd)). 

Conductance plateau features similar to those found 
in Fig.[T]^d) are also observed in other 30-nm-wide GNR 
devices for channel lengths L up to 1.7 /xm. Figs.[2](a) 
and (b) show the G{Vg) measured for two GNRs with 
L = 900 nm and 1.7 /im, respectively. At low temper- 
atures, both GNRs exhibit equally-spaced conductance 
steps as Vg varies. From Fig.[ljd) and Fig. [2] while the 
spacing AG between conductance plateaus depends on 
the channel length, where AG decreases with increas- 
ing L, the number of plateaus is roughly constant in the 
same gate voltage span (^ 20 V) for these GNRs of the 
same width. It is also important to note that these con- 
ductance plateau structures are highly reproducible un- 
der different thermal cycles, and they generally become 
more well-defined as T decreases, in particular for longer 
channel devices. However, as the temperature further de- 
creases below 10 K, the G{Vg) traces of these GNRs are 
usually overwhelmed by a pronounced fluctuation back- 
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FIG. 2: Conductance measured as a function of gate voltage 
for 30-nm-wide GNR devices with different channel lengths. 
The GNR channel lengths are 900 nm and 1.7/im for (a) and 
(b), respectively. At low temperatures, both devices show 
equally-spaced conductance plateaus with AG ~ 1.7 /xS and 
0.6 /iS for (a) and (b), respectively. The bias voltage is kept 
at 10 mV. 

ground as a function of gate voltage. Unlike the plateau 
structures, the low-T background fluctuations are not re- 
traceable in different thermal cycles, indicating that the 
low-T background fluctuations may be due to the univer- 
sal conductance fluctuation (UCF) phenomena |6j . 

Based on these experimental findings, we attribute the 
conductance plateaus observed in GNRs to the forma- 
tion of ID conduction modes in nanoribbons. There have 
been numerous theoretical studies on the GNR electronic 
structure and its dependence on the ribbon width and the 
cutting angle using 7r-orbital tight-binding models [ITJ HI] 
and the first principle calculations |13[ fT5] . To a first or- 
der approximation, the band structure of GNRs can be 
described by zone folding the graphene band structure. 
In a graphene nanoribbon of width W, the wave vector 
perpendicular to the transport direction has a quantiza- 
tion requirement k±W = nm, giving rise to various ID 
subbands, each with the dispersion relation given by 

E„,m) = ±hvf^^ + {m + a)27r2/VF2, (1) 

where m ~ 0,±1,±2,... is an integer for the subband 
index and k± and fcy are the wave- vector perpendicular 
and parallel to the transport direction, respectively. Here 
< |q;| < 0.5 depends on the crystallographic orientation 
of the GNR, and yields a bandgap Agnr = 2A£^|a|, 
where AE = hvfTr/W is the energy separation between 
the ID subbands. We note that in GNRs, the wavefunc- 
tion has to vanish at the ribbon edges, in contrast to the 
periodic boundary condition in the case of carbon nan- 
otubes (CNTs). This leads to a different quantization 
requirement in CNTs, namely k±WcNT — ^nm, and as 
a result twice as large energy separation in CNTs for 
the same circumference length W. In addition, the ID 
subbands of GNRs are singly degenerate as opposed to 
those of CNTs [TO], where the orbital degeneracy asso- 
ciated with the K-K' bands is lifted with the splitting 
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determined by the GNR chirality. 

For a more realistic model of the GNR band struc- 
ture, we adopt the tight-binding method with the next- 
neighbor hopping tc — 2.7 eV. The Fermi velocity is 
related to the hopping integral t according to Vf — 
3atc/2h ~ lO^m/s, where a « 0.246 nm is a lattice 
constant of graphene. At the ribbon edge, we assume 
enhanced hopping between the neighboring C atoms by 
12% to account for the C-C bond length contraction [T7] . 
We find that in armchair GNRs (unrolled zigzag CNTs), 
a « 0.27, 0.4, and 0.066 in the families N = 3p, 
N = 3p-|-l, = 3p-|-2 respectively, where p is an 
integer and N is the number of dimer lines across the 
ribbon width. As a result, the K-K' orbital splitting in 
families N = 3p+l and A = 3p + 2 is about 0.2A£; and 
0.13A£^, respectively. On the other hand, for the N — 3p 
family, the subband energies are separated by 0.46Ai? (or 
0.54A£'), such that the resulting subbands are roughly 
equally spaced with half the quantization energy AE. 

Using the Landauer approach, the device conductance 
at a finite temperature can be expressed by 

a.^^^frmi^'^y^. (2) 

where Ti{E) is the transmission probability of carriers 
in each subband at energy £', and the factor 2 comes 
from the spin degeneracy. At zero temperature, Eq. ([2]) 
can be reduced to G = 2e^ /h^^ti where U — Ti{Ef) 
and the summation includes all the ID modes below the 
Fermi energy. In our simulation, we assume an energy- 
independent transmission probability t for all subbands. 

Now we compare simulations of the GNRs with the 
experimental results. For a 30-nm-wide GNR as shown 
in Fig.[l] we estimate the quantization energy AE to be 
60 meV. From simulations based on Eq. ([2| , we find that 
an energy spacing of at least 6kT is required in order 
to observe conductance quantization features at a given 
temperature T. Since the conductance steps are clearly 
visible in Fig.[l]at 80 K, corresponding to a minimum en- 
ergy separation of 40meV, the energy splitting at the 
K-K' bands must be smaller than 20meV. These facts 
can be accounted for, within the tight-binding model, by 
adopting a (3p-|-l) family armchair ribbon to model the 
electronic structure of the GNR in Fig.[T| Fig.|3|a) shows 
the calculated density of states of an armchair GNR with 
A^ = 253, corresponding to M^=31nm. The quantization 
energy in this ribbon is AE « 58 meV and the separation 
between the nearly degenerate pairs of bands is about 
AE' ~ 0.8AE w 46 meV. Indeed, we find an excellent 
agreement between the measurement and the simulation 
at 80 K as shown in Fig.jsj^b). In Fig.[3](b), the step height 
is Ae^ /ht because the steps associated with individual oc- 
cupation of nearly degenerate bands, which are separated 
by only 11 meV, cannot be resolved at this temperature. 
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FIG. 3: (a) The density of states of a graphene nanorib- 
bon. (b) Calculated normalized conductance G/{2te^ /h) as 
a function of carrier density for graphene nanoribbon devices 
of VK = 31 nm, using armchair GNR with A'' = 253. The inset 
illustrates schematically the field lines and the equipotential 
surfaces of the GNR, demonstrating a significant discrepancy 
from the case of a parallel-plate capacitor. 



In order to further resolve these nearly degenerate sub- 
bands, the measurements should be performed below 10 
K, the regime where universal conductance fluctuation 
dominate the transport properties. From Fig.[3]^b), we 
calcuate the transmission probability t of this GNR to 
be 0.016. We note that while the transmission proba- 
bility may be a function of the length, energy and the 
contact to the sample, the experimental results Fig.|3]^b) 
can be well-described by a constant t at low energies. On 
the other hand, the discrepancy between the simulated 
and measured conductance of the fourth plateau around 
Vg > \7V may be associated with the enhanced trans- 
mission probability at high energies. 

The carrier density in the GNR is related to the gate 
voltage as p = Cg{Vg — Voirac), where Cg is a gate ca- 
pacitance and l^Dirac IS a gate voltage corresponding to 
the Dirac point. As shown by Fig.|3jb), we observe four 
conductance plateaus over a Vg span of 20 V, correspond- 
ing to a shift of the Fermi level by E^ « 260 meV and a 
change in carrier density of 8.4 x 10^^ cm~^. This is com- 
parable to the carrier density in 2D graphene at the same 
Fermi energy Ej given by p2D = Ef/nvf. The depen- 
dence of p on Vg leads to an effective gate capacitance 
of Cgnr ~ 4.2 X lO^^cm^^V^^, which is significantly 
larger than the value expected for 2D graphene, given 
by C2D — 7.2 X 10^°cm~2V~^ for the gate dielectric of 
Si02 with thickness h — 300 nm. This difference is due 
to the fact that when the aspect ratio of the GNR width 
to the oxide thickness (w/h) becomes much smaller than 
one, the field lines between the GNR and the gate devi- 
ates from that of a parallel-plate capacitor (see inset of 
Fig.|3jb)), resulting in an increasing Cgnr/C'2d ratio as 
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FIG. 4: (a) Conductance of a GNR device measured as a 
function of gate voltage at T=15 K for two drain biases, 50 mV 
and 100 mV. The GNR device has a channel length and width 
of 1.7 /im and 30 nm, respectively, (b) Schematic energy band 
diagram of a GNR with a drain bias voltage applied. The 
conductance is determined by the number of ID subbands 
accessible in the bias window. For a given gate voltage, the 
number of ID modes contributing to the transport can be 
increased at a sufficiently large drain bias. 



(w/h) decreases. We have numerically calculated the ca- 
pacitance of the GNR of infinite length L in the back gate 
configuration where the Si02 dielectrics and the vacuum 
below and above the ribbon are explicitly considered. 
For /i/w=300/30, we find Cgnr/Cq ^ 10[T8], which is 
within a factor of 2 of the estimated values based on 
experiments [12] . 

We also note that, since the device conductance at 
room temperature is comparable to its low temperature 
value, the scattering mechanism is not temperature sen- 
sitive in these GNRs, and therefore, the dominate scat- 
tering process is mostly due to the impurity scattering 
within the channel or the edge disorders. This results 
in a transmission probability that decreases with the in- 
creasing channel length. It also implies that these ID 
transport channels are not ideally ballistic. The fact that 
similar conductance quantization is observed for both 
electrons and holes reflects the symmetry in the electron 
and hole band structures as well as in their transmission 
probability. 

In addition to the Fermi level modulation by the gate 
control, the number of ID conduction modes that con- 
tribute to the transport in GNRs also depends on the 
bias voltage between the source and drain contacts, yield- 
ing a Vrf-dependent device conductance in GNRs. [TT] 
Fig.|4|^a) shows measured conductance of a GNR device 
{W = 30 nm and L = 1.7 fj,m) as a function of Vg for 
two bias voltages Vd of 50 mV and 100 mV at 15 K, both 
exhibiting conductance plateaus with the same quantiza- 
tion spacing AG as Vg varies. However, in Fig.|4]^a), it is 
noted that the conductance curve G{Vg) at Vd — 100 mV 
is always higher than that of Vd = 50 mV by roughly one 
unit of quantized conductance AG for all gate voltages. 
In comparison, at Vd = 10 mV, this GNR possesses a 
nearly identical G{Vg) curve to that at Vd = 50 mV. This 
dependence of conductance on the drain bias can be un- 



derstood through the schematic shown in Fig.|4]^b). As 
the source-drain bias windows becomes larger than the 
energy spacing of the ID subbands, which is ^ 67meV 
for a 30-nm-wide GNR, the number of ID conduction 
modes is increased by one regardless of the gate voltage, 
resulting in a higher conductance level by AG. 

In conclusion, we have fabricated graphene nanorib- 
bon devices and reported the first experimental observa- 
tion of conductance quantization phenomena in GNRs 
at temperatures as high as 80 K and channel lengths 
up to 1.7/im. We have performed temperature and 
bias dependence studies, and shown that the energy 
spacing between ID subbands of 30-nm-wide GNRs is 
around 50 meV. The experimental results are in excellent 
agreements with theoretical calculations within the tight- 
binging approximation. The experimental findings here 
provide an important step towards developing graphene- 
based quantum devices. 

The authors thank Bruce Ek for expert technical as- 
sistance. 
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The gate capacitance of GNRs can also be approximated 
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by the following semi-analytical approach. Assuming a 
homogenous dielectric medium e below and above the 
GNR, the capacitance due to the back gate can be ex- 

pressed by ^ = f arctan + In (l + ^) , where 
Co = ee/h. The solution assumes a constant charge den- 
sity on the GNR. The gate capacitance of the device is 
then approximated by Cgnr ^ {Cq" + Cg''~'^)/2, yielding 
an enhancement factor of 9. 
[19] The actual gate capacitance of the GNR is expected to 
be smaller than the numerical simulations because of (i) 
the non-ideal dielectrics used and (ii) the screening from 
source and drain electrodes. 



